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1. Introduction 


The data throughput that a wireless sensor network (WSN) can guarantee is influenced by 
a plethora of concurrent causes. Among those, limited connectivity and medium access 
control (MAC) failures are major issues that should be carefully considered. The aim of this 
chapter is to provide the reader with a neat and general mathematical framework for the an- 
alytical computation of key performance metrics of WSNs. The focus is on connectivity and 
MAC issues. Quantitative answers to such questions as the following will be given: how well 
is the network -or a subset of it- connected? What is the rate at which sensors are able to 
transmit their data to sink(s)? What is the overall throughput of a sensor network deployed 
on a specific domain? 

We consider a multi-sink WSN where sensor and sink nodes are both randomly deployed on a 
finite or infinite domain. Sensors are in charge of sampling the surrounding environment and 
send their data to one of the sinks, possibly the one providing the best signal strength. The 
computation requires some basic assumptions that hold throughout the chapter: two nodes 
are considered connected if the path loss (including both a deterministic distance-dependent 
component and a random fluctuation) is above a fixed threshold; all nodes employ the same 
transmission power; sinks have an ideal connection to an infrastructured processing center. 
We first address connectivity issues by considering single-hop networks with nodes deployed 
on the infinite plane, then, after discussing the role of border effects and providing a mathe- 
matical means to deal with them, we consider networks on finite regions of square shape. The 
probabilities that a randomly chosen sensor is connected to one of the sinks, that all sensors 
-or some percentage of them- are connected, are computed. The connectivity model is then 
generalized to handle the case of rectangular deployment regions as well as inhomogeneous 
nodes densities. However, signal strength based connectivity is not exhaustive for real-life 
applications where failures may occur due to packet collisions, even in perfect channel condi- 
tions. For this reason, we also present a rigorous approach for modeling the MAC layer under 
a carrier-sense multiple access with collision avoidance (CSMA/CA) protocol when several 
sensor nodes compete for accessing the same channel at the same time. In particular, the anal- 
ysis is carried out in the specific case of IEEE 802.15.4 MAC algorithm under both Beacon- and 
Non Beacon-Enabled operation modes. By looking at a single sink scenario with a number of 
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sensors, the practical outcome is the probability of successful packet reception by the sink, 
used to derive the throughput from sensors to sink. 

Finally, going back to a multi-sink scenario, we now have the means for computing the prob- 
abilities that a sensor is connected to an arbitrary sink and that it succeeds in transmitting 
its packet. Therefore, by integrating the two building blocks mentioned before, we end up 
with an analytical tool for studying the performance of multi-sink WSNs, where MAC and 
connectivity issues are taken into account. Network performance is synthesized by introduc- 
ing the concept of Area Throughput, that is, the number of samples per unit of time success- 
fully delivered by the sensors to the infrastructure. Numerical results are given for the case 
of IEEE 802.15.4 MAC protocol. The model is also applicable to WSNs employing any MAC 
protocol. 

The chapter is organized as follows. In Section 2 the application scenario is described and 
some related works are presented. Section 3 introduces the link and connectivity models used. 
In Sections 4 and 5 connectivity results are derived for the case of unbounded and bounded 
networks, respectively. Section 6 is devoted to the MAC model and finally Section 7 reports 
throughput results. 


2. Application Scenario 


A multi-sink WSN is considered where data collection from the environment is performed 
by sampling some physical entities and sending them to some external user. The reference 
application is spatial/temporal process estimation Verdone et al. (2008) and the environment 
is observed through queries/response mechanisms: queries are periodically generated by the 
sinks, and sensor nodes respond by sampling and sending data. Through a simple polling 
model, sinks periodically issue queries, causing all sensors perform sensing and communi- 
cating their measurement results back to the sinks they are associated with. The user, by col- 
lecting samples taken from different locations, and observing their temporal variations, can 
estimate the realisation of the observed process. Good estimates require sufficient data taken 
from the environment. Often, the data must be sampled from a specific portion of space, even 
if the sensor nodes are distributed over a larger area. Therefore, only a location-driven sub- 
set of sensor nodes must respond to queries. The aim of the query/response mechanism is 
then to acquire the largest possible number of samples from the area. Since the acquisition 
of samples from the target area is the main issue for the application scenario considered, a 
new metric for studying the behavior of the WSN, namely the Area Throughput, denoting the 
amount of samples per unit of time successfully transmitted to the final user originating from 
the target area, is defined. As expected, area throughput is larger if the density of sensor 
nodes is larger; on the other hand, if a contention-based MAC protocol is used, the density 
of nodes significantly affects the ability of the protocol to avoid packet collisions (i.e., simul- 
taneous transmissions from separate sensors toward the same sink). In fact, if the number of 
sensor nodes per cluster is very large, collisions and backoff procedures can make data trans- 
mission impossible under time-constrained conditions, and samples taken from sensors do 
not reach the sinks and, consequently, the final user. Therefore, the optimization of the area 
throughput requires proper dimensioning of the density of sensors, in a framework model 
where both MAC and connectivity issues are considered. Although our model could be ap- 
plied to any MAC protocol, we particularly refer to CSMA-based protocols, and specifically 
to IEEE 802.15.4 air interface. In this case, sinks act as PAN coordinators periodically trans- 
mitting queries to sensors and waiting for replies. According to the standard, the different 
personal area network (PAN) coordinators, and therefore the PANSs, use different frequency 
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channels. Therefore no collisions may occur between nodes belonging to different PANss; 
however, nodes belonging to the same PANs compete when trying to transmit their packets 
to the sink. An infinite area where sensors and sinks are uniformly distributed at random, is 
considered. Then, a specific portion of space, of finite size and given shape (without loss of 
generality, we consider a square or a rectangle), is considered as target area (see Figure 1). 


X sensor 
e. sink 


Fig. 1. The Reference Scenario considered. 


We assume that sensors and sinks are distributed over the bi-dimensional plane with densities 
ps and po, respectively, with the latter much smaller than the former. Denoting with A the area 
of the target domain and by k the number of sensor nodes in A, k is Poisson distributed with 
mean k = os - A and p.d.f. 


Be i (1) 


We also let I = pọ - A be the average number of sinks in A. 

The frequency of the queries transmitted by the sinks is denoted as f; = 1/T,. Each sensor 
takes, upon reception of a query, one sample of a given phenomenon and forwards it through 
a direct link to the sink. Once transmission is performed, it switches to an idle state until 
reception of the next query. We denote the interval between two successive queries as round. 
The amount of samples available from the sensors deployed in the area, per unit of time, is 
denoted as Available Area Throughput. In this Chapter we determine how the area throughput 
depends on the available area throughput for different scenarios and system parameters. 


2.1 Related Works 

Many works in the literature devoted their attention to connectivity in WSNs or to the ana- 
lytical study of carrier-sense multiple access (CSMA)-based MAC protocols. However, very 
few papers jointly consider the two issues under a mathematical approach. Some analysis of 
the two aspects are performed through simulations: as examples, Stuedi et al. (2005) related 
to ad hoc networks, and Buratti & Verdone (2006), to WSN. Many papers based on random 
graph theory, continuum percolation and geometric probability Bollobàs (2001); Meester & 
Roy (1996); Penrose (1993; 1999); Penrose & Pistztora (1996) addressed connectivity issues of 
networks. In particular, wireless ad hoc and sensor networks have recently attracted a grow- 
ing attention Bettstetter (2002); Bettstetter & Zangl (2002); Pishro-Nik et al. (2004); Salbaroli & 
Zanella (2006); Santi & Blough (2003); Vincze et al. (2007). A great insight on connectivity of 
ad hoc wireless networks is provided in Bettstetter (2002); Bettstetter & Zangl (2002); Santi & 
Blough (2003). Nonetheless, the authors do not account for random channel fluctuations and 
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do not explicitly discuss the presence of one or more fusion centers (sinks) in the given re- 
gion. Connectivity-related issues of WSNs are addressed in Salbaroli & Zanella (2006); Vincze 
et al. (2007). In Salbaroli & Zanella (2006), while considering channel randomness, the authors 
restrict the analysis to a single-sink scenario. Although single-sink scenarios have attracted 
more attention so far, multi-sink networks have been increasingly considered in the very re- 
cent time. As an example, Vincze et al. (2007) addresses the problem of deploying multiple 
sinks in a multi-hop limited WSN. However, the work presents a deterministic approach to 
distribute the sinks on a given region, rather than considering a more general uniform random 
deployment. Furthermore, since the finiteness of deployment region plays a not secondary 
role on connectivity, those models based on bounded domains turn out to be of more practical 
use. 

Concerning the analytical study of CSMA-based MAC protocols, in Takagi & Kleinrock (1985) 
the throughput for a finite population when a persistent CSMA protocol is used, is evaluated. 
An analytical model of the IEEE 802.11 CSMA-based MAC protocol, is presented by Bianchi 
in Bianchi (2000). In these works no physical layer or channel model characteristics are ac- 
counted for. Capture effects with CSMA in Rayleigh channels are considered in Zdunek et al. 
(1989), whereas Kim & Lee (1999) addresses CSMA/CA protocols. However, no connectivity 
issues are considered in these papers: the transmitting terminals are assumed to be connected 
to the destination node. In Siripongwutikorn (2006) the per-node saturated throughput of an 
IEEE 802.11b multi-hop ad hoc network with a uniform transmission range, is evaluated un- 
der simplified conditions from the viewpoint of channel fluctuations and number of nodes. 
Also, some studies have tried to describe analytically the behavior of the 802.15.4 MAC pro- 
tocol. Few works devoted their attention to non beacon-enabled mode (see, e.g. Kim et al. 
(2006)); most of the analytical models are related to beacon-enabled networks Misic et al. (2004; 
2005; 2006); Park et al. (2005); Pollin et al. (2008). Some of these fail to match simulation results 
(see, e.g. Pollin et al. (2008)), whereas slightly more accurate models are proposed in Park et al. 
(2005) and Chen et al. (2007), where, however, the sensing states are not correctly captured by 
the Markov chain. In conclusion, the most relevant difference between the previously cited 
models and the one developed in Buratti & Verdone (2009) and Buratti (2009) and used here, 
is that the latter precisely captures the algorithm defined by the standard, while considering a 
typical WSN scenario. In our scenario nodes only have one packet to transmit to the sink (i.e., 
when they receive the query and have to transmit data before the reception of the subsequent 
query). Therefore, the number of nodes competing for channel at a given time is unknown 
and not constant (as it is in the above cited works) but it decreases with time, since successful 
nodes go to sleep till next query. 

Finally, to the best of the Authors knowledge, no one has so far introduced any 
connectivity /MAC model for WSNs while jointly considering the following aspects: pres- 
ence of both sensors and multiple sinks, random deployment of nodes, bounded scenarios, 
channel fluctuations, realistic MAC protocol in non-saturation condition. 


3. Link and Connectivity Models 


Many works in the WSN scientific literature assume deterministic distance- dependent and 
threshold-based packet capture models. This means that all nodes within a circle centered at 
the transmitter can receive a packet sent by the transmitting one Bettstetter (2002); Bettstet- 
ter & Zangl (2002); Santi & Blough (2003). While the threshold-based capture model, which 
assumes that a packet is captured if the signal-to-noise ratio (in the absence of interference) 
is above a given threshold, is a good approximation of real capture effects, the deterministic 
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channel model does not represent realistic situations in most cases. The use of realistic channel 
models is therefore of primary importance in wireless systems. 

In this chapter, a narrow-band channel, accounting for the power loss due to propagation 
effects including a distance-dependent path loss and random channel fluctuations, is consid- 
ered. 

Specifically, the power loss in decibel scale at distance d is expressed in the following form 


L(d) = kg + kı Ind + s, (2) 


where kg and kı are constants, s is a Gaussian r.v. with zero mean, variance a2, which rep- 
resents the channel fluctuations. This channel model was also adopted by Orriss and Barton 
Orriss & Barton (2003) and other Authors Miorandi & Altman (2005). In Verdone et al. (2008) 
experimental measurement results, performed with 802.15.4 devices at 2.4 [GHz] Industrial 
Scientific Medical (ISM) band, deployed in different environments (grass, asphalt, indoor, etc), 
are shown. It is found for the received power in logarithmic scale that in general a Gaussian 
model can approximate the measurement variation fairly well, with different values of the 
standard deviation. By suitably setting kı, it is possible to accommodate an inverse square 
law relationship between power and distance (k4 = 8.69), or an inverse fourth-power law 
(kı = 17.37), as examples. 

For what concerns the link model, a radio link between two nodes is said to exist, which means 
that the two nodes are connected or audible to each other !, if L < Li, where Ly, represents the 
maximum loss tolerable by the communication system. The threshold Li, depends on the 
transmit power and the receiver sensitivity. 

By solving (2) for the distance d with L = Lip, we can define the transmission range 


Lg, —-ko—5 


TR=e ^ , (3) 


as the maximum distance between two nodes at which communication can still take place. 
Such range defines the connectivity region of the sensor. Note that by adopting independent 
r.V.'s s for separate links, we have different values of TR for different sinks, given a generic 
sensor. In other words, unlike many papers dealing with connectivity issues in the literature 
Bettstetter (2002); Bettstetter & Zangl (2002); Santi & Blough (2003), we do not use circles to 
predict sensor connectivity. However, by setting 7 = 0, we neglect the channel fluctuations 
and may still define an ideal transmission range, as a reference, as 


Len -ko 


ZR | N (4) 
Finally, we can define a connection function between any node pair whose distance is d as 


o 


g(d) = Prob {L(d) < Ly} =1- Lerfe (= — WD (5) 


3.1 Connectivity properties in Poisson fields 

Connectivity theory studies networks formed by large numbers of nodes distributed according 
to some statistics over a limited or unlimited region of R, with d=1,2,3, and aims at describing 
the potential set of links that can connect nodes to each other, subject to some constraints from 
the physical viewpoint (power budget, or radio resource limitations). 


1 Jink's reciprocity is assumed. 
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It is widely accepted that, a WSN is fully-connected in case any sensor node is able to reach at 
least one sink node, either directly or through other sensor nodes Verdone et al. (2008) (not 
necessarily requiring any node to be reached by any other node). 

Let us consider a stationary Poisson Point Process (PPP) ® = {x1,x2,...} having intensity 
p, with x; = (xi, yi), i = 1,2,... being a random point in R*. ® may also be regarded as 
a random measure on the Borel sets in IR*: taken any Q C IR? having area Wo, ®(O) is a 
Poisson r.v. which counts the number of points of o that lie in the set O, whose first order 
moment is 


E(&(Q)) = pu (O) = p | dx = pWa, © 


where v4 (Q) is the Lebesgue measure of O. Now suppose we want to count only those points 
in © which are connected to an arbitrary node xg: this implies a thinning procedure on ® 
such that each point is retained with probability C(||xo — x;||) and discarded with probability 
1 —- C(||xo —»x|]), i = 1,2,..., where C(x) is a non-negative measurable function such that 
0 € C(x) € 1. By so doing, the new inhomogeneous process ®' is obtained. 

By recalling the Campbell Theorem for point processes Gardner (1989) that we report for later 


E( E fo) =o [fait O 


xEQ 
for any non-negative measurable function f, we have for ®’ 


p = E($'(Q)) =E (E C(|Xo — D = e]. C(||xo — x||)dx. (8) 
xEQ 

In particular, when the channel model of eq. (2) is used (i.e., C(x) = g(x)), the mean number 

of nodes audible within a range of distances r; and r, to a generic node (r > r1), is denoted as 

Hr," and can be written as Orriss & Barton (2003); Orriss et al. (1999) 


Un r = np|Y (ay bir) — F(a, biri), (9) 
where p is the initial nodes’ density and 
V(a,b;r) = r20(a,—b,Inr) 
2a, + p. (10) 
— ad * (n — bı Inr+2/b,), 


and a4 = (Lg, — ko) /o, by = ky /o and ®(x) = [3,07 2n)e-" du. 


4. Connectivity in Unbounded Networks 


Since the channel model described by eq. (2) is used, the number of audible sinks within a 
range of distances rı and r from a generic sensor node (r > 71), n,,;, is Poisson distributed 
with mean j/;, ;, given by eq. (9) by simply substituting p with pg. Then by letting rı = 0 and 
r — co, we obtain 

Hoo = Ttpo exp[(2(Ltn — ko) /k1) + (207 /Kr)]. (11) 
Equation (11) represents the mean value of the total number, no œ, of audible sinks for a 
generic sensor, obtained considering an infinite plane Orriss & Barton (2003). 
Its non-isolation probability is simply the probability that the number of audible sinks is 
greater than zero 
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Jo =1—e P", (12) 


5. Connectivity in Bounded Networks 


When moving to networks of nodes located in bounded domains, two important changes 
happen. First, even with pọ unchanged, the number of sinks that are audible from a generic 
sensor will be lower due to geometric constraints (a finite area contains (on average) a lower 
number of audible sinks than an infinite plane). Second, the mean number of audible sinks 
will depend on the position (x, y) in which the sensor node is located in the region that we 
consider. The reason for this is that sensors which are at a distance d from the border, with 
d ~ TR;, have smaller connectivity regions and thus the average number of audible sinks 
is smaller. These effects, known in literature as border effects Bettstetter & Zang] (2002), are 
accounted for in our model. 
The result (9) can be easily adjusted to show that the number of audible sinks within a sector of 
an annulus having radii rı and r and subtending an angle 26, is once again Poisson distributed 
with mean 

Urro = Opo[¥ (a1, bi; r) — Y(2,b;11)], (13) 


0 < 0 < zt. If the annulus extends from r to r + dr, and 0 = 0(r), this mean value becomes 


óY (a4, b1;r 
iera orp tei oy, 0<89< r. (14) 


Consider now a polar coordinate system whose origin coincides with a sensor node. As a 
consequence of (14), if a region is located within the two radii r; and rz and its points at a 
distance r from the origin are defined by a 0(r) law (see Fabbri & Verdone (2008), Fig. 1), 
then the number of audible sinks in such a region is again Poisson distributed with mean 


Hrs ze 2 (y Vega dy, that is, from (10) and after some algebra, 
T2 
Urro) = [ 20(r)oor® (a; — bı Inr)dr. (15) 


5.1 Square Regions 

Now consider a square SA of side L meters and area A = L?, sensors and sinks uniformly 
distributed on it with densities ps and po, respectively. Equation (15) is suitable for expressing 
the mean number of audible sinks from an arbitrary point (x,y) of SA, provided that such 
point is considered as a new origin and that the boundary of SA is expressed with respect to 
the new origin as a function of ri, r2 and @(r). In order to apply equation (15) to this scenario 
and obtain the mean number, u(x,y), of audible sinks from the point (x,y), it is needed to 
set the origin of a reference system in (x, y), partition SA in eight subregions (5,1 ...5,8) by 
means of circles whose centers lie in (x, y) (see Fabbri & Verdone (2008), Fig. 2). Thank to the 
properties of Poisson r.v.’s, the contribution of each region can be summed and we obtain an 
exact expression for 


u(x, y) a 20;(r) - po -7- (a — by lnr)dr, (16) 


which is the mean number of sinks in SA that are audible from (x, y), where r, j, r? ;, 0;(r) are 
reported in Fabbri & Verdone (2008), Tables 1-2. 
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If we assume a single-hop network, a sensor potentially located in (x, y) is isolated (i.e., there 
are no audible sinks from its position) with probability p(x, y) = e^! *) and it is non isolated 
with probability 


g(x,y) =1- e M Gw) (17) 


Owing to the assumption that sensor nodes are uniformly and randomly distributed in SA, if 
we now want to compute the probability that a randomly chosen sensor node is not isolated, 
we need to take the average q(x,y) on SA. In fact, the probability that a randomly chosen 
sensor node is not isolated (which is an ensemble measure) and the average non-isolation 
probability over a single realization coincides due to the ergodicity of stationary Poisson pro- 
cesses (see Stoyan et al. (1995), page 104). This was also verified by simulation. 

Recalling that we have considered the lower half of the first quadrant, which is one eighth of 
the totality, we have 


_ 8 [L/2 rx 
qe f 1 q(x, y)dydx. (18) 


5.2 Rectangular Regions 

We now consider a rectangular domain C of sides S4 and 55, S4 > S5, area W = 5, S5, with 
sensors and sinks uniformly distributed on it with densities ps and po, respectively. We aim at 
computing the mean number of audible sinks from a fixed position (x, y) which are contained 
in C. Since we are dealing with a rectangular domain whose points have to be expressed in 
polar coordinates in order to apply (15), such a domain has to be properly partitioned into a 
set of subregions, to be defined in terms of r1, r2, and 0. Moreover, unlike the case of square 
domain, the nature of the partition depends on the position (x, y) considered. In particular, if 
we restrict the analysis to the upper-right quart, we can identify 4 different cases depending 
on whether (x, y) belongs to A1, A», A3 or A4 (see Figure 2). Let us denote as case i the event 
(x,y) € Aj, for i — 1,2,3,4. In each of the latter cases, the domain is differently partitioned 
into 8 subregions that are sectors of annuli. What changes from one case to another is the 
definition of each subregion. As an example, the subregion having r in the range [0, 51 /2 — y| 
lies completely in C only when (x, y) € A»; otherwise it partially exceeds the borders of C. 
Thus, the corresponding angle 0(r) is zt in case 2 and some function of r in the other cases. The 
following tables define A1-A4 and the values of r and 0 in each subregion for case i = 1,2,3,4, 


respectively. In the following, we denote by PDA the range of r of the jth subregion 


when in case i, and by fa ) (r) the corresponding angle. 
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(0, 91/2) 


(0, —54/2) 


(Se, 54/2) 


C (S2, —51/2) 


Fig. 2. Geometric partitioning of the rectangular region. 


(xy)|(91/2€x€ $5, 0< y< x- 51/2) 


(x,y) | {52/2 < x < S2, x 51/2785 < y € 54/2] 


(xy) | {52/2 <x € $5, max($4/2-x,x — 81/2) < y € $0/2—85 +x} 


(x,y) | {52/2 < x < $4/2, 0€y € $/2- x) 


(41) 


Range: 74 


(Ay) 
<r< r3 


O<r<Sy-x 
Sy -—x<r<S1/2-y 
$1/2-y€r« (S5 —x)2 + ($4/2- y? 
V (Sq 3)? + (S4/2- y < r < 84/2 y 


(S2 - 2 + (51/2 +y)? 


Sy/2+y<r< 
(Sy —x)? (55/2 y? Sr<x 


qc x2 + ($41/2- y? 


y2 +(S1/2- y)? <r < y2 t (54/24 y)? 


Region Range: "Un <r< 2) 
1 0zr«54/2-y 
2 $1/2-y €&r«$5-x 
3 S2 =x <r < y (S3 - x) + (51/2 — y)? 
4 VS 32 + (5/27 3 <r <x 
5 rare + (S1/2- y? 


6 V3? + (S1/2-y)2 <r < S1/2+y 


7 S1/2+y <r < ($5 x)? + (81/2 + y? 


8 V (52 33 € (4/24? <r < yx2+( 
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Region Range: "Ma Er« 763) e(43) (r) 
1 0<r<S-x mT 
2 Sy -x <r<S1/2-y Z + arcsin S2 
3 Sy/2-y<r< (Sq — x)2 + ($4/2 — y)2 H + arcsin LOY — arccos S23 
4 ($5 — x)2 + (51/2— y?zr«x 5 | 1 arcsin Spy — arccos Buc 
5 yere NE + (S4/2- y)? 5 — arccos Bey + 1 (arcsin 81 foc — arccos — 


6 V3 + (S1/2-y)2 <r < S1/2+y 5 


č $5—x 
— i (arccos x + arccos 2 


$5-x 


7 Sy/2+y <r < J (S5 x)? + (4/24 y? arcsin Sit + 4 (arccos 7— +arccos & 
8 V2 - 32 (55/24? <r < V2 (51/24 y? i arcsin X — arccos -— 


Region Range: Ko <r< fia) g(Aa) (r) 
1 O<r<S)-x ps 
1 _ $9—X 
2 Sg-x<r<x 5 + arcsin —; 
3 x€r«$1/2-y F + arcsin 52—* arcos $ 
4 5$1/2-y&r« (Sq — x)2 + ($4/2- y? 5 + arcsin S1/2Y — arccos € — arccos pct 
$4/2 $4/2— $5—23 
5 V(52 33 + (8,/2- 3 <r < V2 e (42-9)? H — arccos 1/273 41 (arcsin 1/5 Y arcos 2. =) 
6 V3 + (S1/2- y)? <r < S1/2+y F = § (arccos 23 -arccos 8) 
7 Sy/2+y <r < J(S5 x) + (61/2 +y)2 $ — arcos 51/2 = E (arccos # + arccos meu) 
8 V2 - 2 (55/24? <r < V3 (51/24? 1 arcsin S1+Y arcos ¥ 


Note that when S; = S; the partitioning scheme degenerates to the one for square regions. 
Now, starting from (15) and owing to the linearity of Poisson independent r.v.’s, the mean 
number of sinks that are audible from (x, y), with (x, y) € Aj, may be computed as 


Ai B. ogg 
uA) (x,y) = x 
jae 


) 
2098) - go: r- (a — bj Inr)dr, 


(19) 


for i = 1,2,3,4 and with a, = (Lg, — ko) /c, by = ky/e and ®(x) = [*,(1/V2n)e-" "?du. 
Owing to the Poisson distribution of the number of audible sinks, the probability that the 
position (x, y), with (x, y) € Aj, is isolated (i.e., no sink is heard) is simply 


p^? (x, y) = e ^? Gan, 


(20) 


while the probability that the position (x, y), with (x, y) € Aj, is not isolated is 


q 9 (x, y) = 1— p (xy) = 1— ec Gn. 


(21) 


Now, the mean number of sinks that are audible from (x,y), with (x,y) € {A1 U A? U A3 U 


Ag}, is 
wan Gey) (x,y) € Ay 
_) Be (xy) uy) € Ao 
x,y) = (22) 
MOM) TY Mey) yE As 
ud (x,y)  ,Qy)eA 
Equally, the isolation and non-isolation probabilities may be computed as 
pAd(xyy) =e) (x,y) EAr 
(A2) = poe? (xy) cA 
p(x,y) = j A (xy) ^ (A3 (x4 (Qs) 5 (23) 
pAsixy) =e HH (x,y) € Ag 
p^) (x y) = e I9 (xy) P (x,y) € A4 
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and 
g(x,y) =1-e HMM)  ,(xy)eA 
(42) (x,y) =1—e y) (xy) EA 
q(x, y) = o y) E 7 (As) (x ( y) 2 (24) 
g3 (x,y) =1= eF VM (x,y) E€ A3 
gaye (x,y) € A, 


respectively. Hence, the average probability of non-isolation over C is 


d 51/2 
4 |F TE xy) = rhah q(x, y)dydx 
i s q ^? (x, y)dyd e dyd 
om So+x (As) dud 51/2 p$,1/2—x (Ay) "T x 
a , + 3" 4 X, X 
È max( ($,/2-xx-8/2) 1 A y) yax $,/2 Jo q ( y) y ) ) 


5.3 Composite Domains 


D 
Fig. 3. Reference scenario for the analysis of composite domains. 


The scenario that we now want to analyze is of the kind of the one depicted in Figure 3. Con- 
sidera D e Pe C' of area W' which is composed of n rectangular sub-domains C’; 


of sides s and s (note that st i) > se holds), area wi), i = 1,2,...,n. We assume the 
sinks are uniformly and fandomiy distributed in C’; with density oo i, i = 1,2,...,n. Instead, 
sensors are uniformly and randomly distributed over the whole domain (i.e., in C^) with den- 
sity ps. As a consequence, sinks are distributed according to a inhomogeneous PPP over C', 
while sensors are distributed according to a homogeneous PPP over C'. 
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Our final goal is to compute the probability that a randomly chosen sensor in C’ is not isolated. 
Now suppose there is a sensor node, S, located in (Sx, Sy) € C', and we want to find the 
probability that it is not isolated. It is clear that the number of sinks that S can hear is not 
limited to the number of sinks contained in C',. Rather, the more its transmission range is 
large compared to the sides of C’;, the more it can benefit from the connectivity offered by the 
sinks located in the other sub-domains (e.g., the adjacent ones). On the contrary, when S is 
not close to one of the borders of C’; and its transmission range is small (i.e., the connectivity 
area of S lies entirely within C',), what happens in C';, Vj # k is totally negligible. We can 
intuitively state that the same happens when po > po, Vj # k, since the other sub-domains 
present too few sinks to provide connectivity to a sensor in C’,. 

Thus, when we are allowed to neglect the interaction between different sub-domains, we can 
simply treat each of them in a separate way. In this way we end up with the n-tuple q = 
(di, dz, ..., jn). The overall approximated non-isolation probability over C’ is obtained as the 
weighted average of q. This case is detailed in Subsection 5.3.1. 

As an alternative, a direct application of (8) with a careful choice of O (i.e., without partition- 
ing) would lead to an exact result. However, the complexity of carrying out the integration 
can sometimes make this approach unfeasible. The details can be found in Subsection 5.3.2. 


5.3.1 Approach 1 
We have q = (91, do, ..., n), with (from (25)) 


PE 
Gi = Ex,y|qi(x,y)] v he ae (x, y)dydx, (26) 


where q;(x, y) is computed on C';, which has sides sË ) and , ) with sË i > sp, peq95./e 
Now, the probability, 7p, that a randomly chosen sensor in C’ is not isolated is simply 


ES i. 
īp = » gi. (27) 


5.3.2 Approach 2 
From (8) and owing once again to the fact that the sum of Poisson independent r.v.'s having 
mean À;, i = 1,2,..., is still Poisson with mean A = Aj + A2 + .. ., we have 


Moro) = Leo fy c (IIx — xol])ax (28) 


i.e., the average number of audible sinks from (xo, yo). 
Equation (28) is very general and takes the interaction between sub-domains into account. 
Now, in order to obtain a result which is analogous to (25), we let 


PM(xo,yo) = e Goin) (29) 


and 

qu (Xo, yo) = 1— pm(xo, yo) = 1 — e Con) (30) 
to end up with the isolation and non-isolation probabilities of the location (xo, yo), respec- 
tively. Then, we simply take the average over the points (xo, yo) such that (xo, yo) € C’ and 
get 


1 
iM = "2 qM (Xo, yo) dxod yo. (31) 
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5.4 Practical Cases With Numerical Results 


5.4.1 Single Rectangle 
Equation (25) can be evaluated numerically once S1, S5, po, Lin, ko, kı, c are known. As an 
example, in Fig. 4 we plot 7 as a function of the ratio y = 55/5. 


e 
[ 


probability 


m 
m 


0.1 : poi Piiii : por prin 
0.001 0.01 0.1 1 


Y 
Fig. 4. 4 as a function of y for different values of Lin, with W = 1 Km, po = 100/W, ko = 40, 
ky = 15,03, 0 = 3. 


As y varies from 1 to 0, the area W remains constant while the domain C gets increas- 
ingly squeezed. The general trend suggests that the smaller is y, the smaller is the level of 
connectivity. This is due to border effects: when S becomes comparable with the transmission 
range, the connectivity area of the sinks is very likely to overstep the domain area, thus re- 
sulting in a decrement in the average number of connected sensors per sink. In particular, we 
expect this to be more appreciable for greater transmission ranges. In fact, from Fig. 4 we can 
observe that for Ly, = 80 dB (TR; ~ 21.54 m), when y ranges from 1 to 0.001 (S2 ranging from 
1000 m to 31.62 m) the loss in connectivity is only (Li = 80 dB; y = 1) — (Lg, = 80 dB; y = 
0.001) ~ 0.04. Instead, for Ly, = 100 dB (TR; ~ 99.96 m) and y ranging as above, the loss in 
connectivity is no less than g(Li, = 100 dB; y = 1) — ĝ(L = 100 dB; y = 0.001) ~ 0.51. 


5.4.2 Composite Domain 
Consider now the non-isolation probability for the composite domain of Figure 3. Assume 


S = 850 m, SU = 400 m, SP = 150 m, s) = 700 m, S) = 400 m, s( = 300 m and the 
densities pọ = 4.E-4, 092 = 3.E-3, 093 = 1.E-3, 00,4 = 6.E-4. 
From (27), the computation of jp is straightforward. In Figure 6 we report jp, 41, 92, 93, da. 
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As for gm, set the origin in D and let (xo, yo) be a generic point in C’. Accounting for the 4 
different zones, the mean number of audible sinks from (xo, yo) is 


4 
MG) = È pox f, Clx- xol Dax (32) 
k=1 k 
Sx ps45M —yo 
= m]. [. C(y/ (x= x0)? + (y — yo)?)dydx 
=X0 Si” —yo 
$0 rs. (pss -yo 5 
+ p02 fi Yo Cy (x — x0)? + (y — yo)?)dydx 
Si —X0 Si —Yo 
54) 4.50) su 


2 +5; —xXo 
teal o [CQ = 20) + Y= yo)?)dydx 
S —X0 Vo 


2 
six psf- 

tma. fi, C Go x) + (y— yo)?)dydx, (33) 
^40 — 40 


while the probabilities of non-isolation of the position (xo, yo) is obtained as 


qM (Xo,yo) = 1 — e "Gon, (34) 


In Figure 5 qu(xo, Yo) is reported. Note that we have qgm(xo, y0) # 0 on the boundaries, a 
fact that confirms that we are not introducing factitious border effects between different sub- 
domains. Note also that equations (32), (33) contain a double integral: this implies a greater 
computational complexity with respect to (19) employed in the Approach 1. On the other 
hand, (32) and (33) are exact (i.e., interactions among sub-domains C"; are not neglected). 
Now, accordingly to (25), the average probability 7,, that a sensor randomly chosen in C’ is 
not isolated is 


sP sP m 
0 


Im = Exp dw Xo. yo)] = n qM (Xo, Vo) dyodxo. (35) 


0 
In Figure 6 we also plot dw as a function of Lin [dB]. It is possible to compare the non-isolation 
probabilities obtained through the two different approaches (bold curves): Approach 2, as 
said, accounts for interactions between sub-domains and thus does not introduce border ef- 
fects that would be fake. This is the reason why we observe gw > Gp (i.e., the WSN performs 
better). Thus jp is a lower bound. 
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Probability 


Fig. 5. dw(xo, yo) on the domain of Figure 3 obtained with L;, = 90 [dB], ko = 40, kı = 13.03, 
g = 3.5. 
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Fig. 6. Non-isolation probabilities referred to the scenario of Figure 3 obtained with kg = 40, 
ky = 13.08, 7 = 3.5. 


6. The IEEE 802.15.4 MAC protocol 


When dealing with contention-based MAC protocols, there exists a certain probability that 
a node does not succeed in accessing the channel or in transmitting its packet correctly (i.e., 
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without collisions). A single-sink scenario, where n 802.15.4 sensors transmit data to the sink 
through a direct link is accounted for, in this Section. We assume all sensor nodes are audible 
to the sink. 

Both, Beacon- and Non Beacon-Enabled modes are considered. We assume that nodes trans- 
mit packets having a size, denoted as z, equal to D - 10 bytes, where D is an integer parameter. 
We also assume that the size of the query packet is equal to 60 bytes.We denote as T the time 
needed for transmitting 10 bytes. Since a bit rate of 250 kbit/sec is used, T = 320psec. 

The Non Beacon-Enabled mode is based on CSMA/CA protocol to access the channel, 
whereas in the Beacon-Enabled case both contention-based and contention-free protocols, are 
implemented. In the latter case a superframe is defined, which starts with a packet denoted 
as Beacon (it coincides with the query packet in our scenario), and divided into two parts: 
inactive and active part. The active part is composed of the Contention Access Period (CAP), 
where a CSMA/CA protocol is used, and the Contention Free Period (CFP), where a max- 
imum number of 7 Guaranteed Time Slots (GTSs) could be allocated to specific nodes (see 
Figure 7, below). The use if GTSs is optional. 

The duration of the whole superframe and of its active part depends on the value of two in- 
teger parameters ranging from 0 to 14, called superframe order, denoted as SO, and beacon 
order, denoted as BO, with BO > SO. In particular, the interval of time between two succes- 
sive Beacons, that is the query interval T; in our scenario, is given by: T; = 16 - 60 - 280. T, 
where T; = 16 usec is the symbol time. Instead, the duration of the active part, denoted as T4, 
is given by: T4 = 16-60- 259 . T,, where 60 - 2°°T; is the slot size. 

The inactive part of the superframe is generally used when tree-based or mesh topologies are 
applied; here, since we are dealing with star topologies, we set SO = BO and T4 = Ty. 

Each GTS must contain the packet to be transmitted and an inter-frame space equal to 40 Ts. 
This is, in fact, the minimum interval of time that must be guaranteed between the reception 
of two subsequent packets. The sink (PAN coordinator, in 802.15.4 jargon) may allocate up 
to seven GTSs; however, a sufficient portion of the CAP must remain for contention-based 
access. The minimum CAP size is 440 T;. By varying packet size D and SO (i.e., the slot 
duration), the number of slots occupied by each GTS and the maximum number of GTSs that 
could be allocated to ensure a CAP larger than 440 T;, will vary as well. As an example, if 
D = 2 and SO = 0, two slots are needed for a GTS, to contain the packet and the inter-frame 
space and a maximum number of 4 GTSs could be allocated. In case SO — 2, instead, each 
GTS will occupy one slot and seven Guaranteed Time Slots (GTSs) could be allocated. We 
denote as Nc s the number of GTSs allocated. 

We assume that in case a node does not succeed in accessing the channel by the end of the 
superframe (in the Beacon-Enabled case) or till reception of the subsequent query (in the Non 
Beacon-Enabled case), the packet will be lost.This implies that by increasing the superframe 
duration the success probability for a node will increase since the node will have more time to 
try to access the channel. Note that in the Beacon-Enabled case, T; may assume only a finite set 
of values (depending on the values of BO); instead, in the Non Beacon-Enabled case T; may 
assume any value. Note that, being (120 + D) - T the maximum delay with which a packet 
can be received by the sink Buratti & Verdone (2009) and having set the query size equal to 
60 bytes, the sink should set T; > (126 + D) - T to make sure all nodes have completed the 
CSMA/CA algorithm. In case lower values of T; are set, a node may receive a new query 
while still trying to access the channel, this resulting in the loss of the old packet. 

We parametrized the behavior of 802.15.4 MAC protocol by means of a function, PyAc(), 
which returns the probability that a sensor node is successful in transmitting its packet when 
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(n — 1) more sensors are trying to do the same. We refer to Buratti & Verdone (2008; 2009) and 
Buratti (2009), Buratti (2010) for derivation and expression of PyAc(n) in Non Beacon- and 
Beacon-Enabled cases, respectively. A finite state transition diagram has been used to model 
sensor nodes states, in both cases Beacon- and Non Beacon-Enabled mode. Here we do not 
report equations for the sake of brevity. In these papers details on formulae are given and also 
a validation of the model against simulation is provided for n < 50 and different values of D. 


6.1 Numerical results 

Some examples of results obtained through the mathematical model developed are shown, 
with the aim of comparing those achieved with the two operation modes (i.e., Beacon- and 
Non Beacon-Enabled). 

In Figures 8(a) PAc(n) as functions of n for the Beacon-Enabled case, for different values of 
SO, with D = 2, is shown. The cases of no GTSs allocated and Ngrs equal to the maximum 
number of GTSs allocable, are considered. As explained above, this maximum number de- 
pends on the values of D and SO. As we can see, PyAc decreases monotonically (for n > 1 
when Ngrs = 0 and for n > Ners when Nets > 0), by increasing n, since the number of 
sensors competing for the channel increases. Once we fix SO, by increasing Ners, Pmac also 
increases, since less nodes have to compete for the channel. Moreover, once Ngrs is fixed, by 
increasing SO, Py Ac also grows, since the CAP size is greater and nodes have a larger amount 
of time to try to access the channel. 

In Figure 8(b) Pmac(n) for different values of D and T;, considering a Non Beacon-Enabled 
network, is shown. As we can see, a decrease of T}, results in a decrement of PyAc, since 
nodes have a smaller amount of time to access the channel. 


Non BE mode 
Query Query Query 
+ >4 > 
Ta Ty 
BE mode 
Beacon/ Beacon/ 
Query Q 
CAP - CFP , Query 


Narts GTSs allocated 


< > 


SD =T, 
Fig. 7. Above part: The IEEE 802.15.4 Non Beacon-Enabled mode. Below part: The IEEE 
802.15.4 Beacon-Enabled mode. 
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Fig. 8. (a): Pmac(n) as a function of n, in the Beacon-Enabled case, for different values of SO 
and Ngrs, having fixed D = 2. (b): Pmac(n) as a function of n, in the Non Beacon-Enabled 
case, for different values of T} and D. 


If we compare the above Figures, we notice that once the superframe duration is fixed, re- 
sults are approximatively the same if no GTSs are allocated, whereas, there is a considerable 
increment of PyAc(n) in the Beacon-Enabled case when GTSs are allocated. Note that the 
cases T, = 15.36 [ms], T; = 30.72 [ms] and T, = 61.44 [ms] correspond to SO = 0, 1 and 2, 
respectively. 


7. Evaluation of the Area Throughput 


The area throughput is mathematically derived through an intermediate step: first the prob- 
ability of successful data transmission by an arbitrary sensor node, when k nodes are present 
in the monitored area, is considered. Then, the overall area throughput is evaluated based on 
this result. 
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7.1 Joint MAC/Connectivity Probability of Success 

Let us consider an arbitrary sensor node that is located in the observed area A at a certain 
time instant. The aim is computing the probability that it can connect to one of the sinks 
deployed in A and successfully transmit its data sample to the infrastructure. Such an event 
is clearly related to connectivity issues (i.e., the sensor must employ an adequate transmitting 
power in order to reach the sink and not be isolated) and to MAC problems (i.e., the number 
of sensors which attempt at connecting to the same sink strongly affects the probability of 
successful transmission). For this reason, we define P; (x,y) as the probability of successful 
transmission conditioned on the overall number, k, of sensors present in the monitored area, 
which also depends on the position (x, y) of the sensor relative to a reference system with 
origin centered in A. This dependence is due to the well-known border effects in connectivity 
Bettstetter (2002). 

In particular, 


Pak(x,y) = EnlPmac(”) - Pcon(x,y)] 
= En[Pmac(n)|- Pcon(x,y)- (36) 


where the impact of connectivity and MAC on the transmission of samples are separated. A 
packet will be successfully received by a sink if the sensor node is connected to at least one 
sink and if no MAC failures occur. The two terms that appear in (36) are now analysed. 
Pcon(x,y) represents the probability that the sensor is not isolated (i.e., it receives a suffi- 
ciently strong signal from at least one sink). This probability decreases as the sensor ap- 
proaches the borders (border effects). Peon for multi-sink single-hop WSNs, in bounded and 
unbounded regions, has been computed in the previous Sections. In particular, for unbounded 
regions, Pcow(x, y) ^ Pcon, that is equal to qo, given by eq. (12). Whereas, when bounded 
regions are considered, Pcow(x, y) is equal to q(x, y) given by eq. (17). 

Specifically, since the position of the sensor is in general unknown, P; (x, y) of (36) can be 
deconditioned as follows: 


Polk = Ex, [Pj (x, y)] 
= Exy(Pcon(x,¥)] -EnlPmac(”)] - (37) 


Ex,|Pcon(x,y)]| is equal to 7 given by, e.g., eq. (25) when a rectangular region is accounted 
for. When, instead border effects are negligible, Ex y|Pcon(x,y)] = Ex,y[Pcon] = Pcon, given 
by eq. (12). 

Given the channel model described in (2) (and following), the average connectivity area of the 
sensor, that is the average area in which the sinks audible to the given sensor are contained, 
can be defined as 


Ag, = ne HW" eh. (38) 

In Fabbri & Verdone (2008) it is also shown that border effects are negligible when Ag, « 0.1A. 
In the following only this case will be accounted for. Thus we have 

Poon (x,y) ~ Pecon 2 1— e ^, (39) 


where uo = PoAc, = IAgz,/ A is the mean number of audible sinks on an infinite plane from 
any position Orriss & Barton (2003), being I = pọ - A the average number of sinks in A. 

Py Ac(n), n > 1, is the probability of successful transmission when n — 1 interfering sensors 
are present introduced in Section 6 for the 802.15.4 MAC case. 
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In general, when CSMA-based MAC protocols are considered, PyjAc(r) is a monotonic de- 
creasing function of the number, n, of sensors which attempt to connect to the same serving 
sink. This number is in general a random variable in the range [0, k]. In fact, note that in (36) 
there is no explicit dependence on k, except for the fact that n < k must hold. Moreover in our 
case we assume 1 < n < k, as there is at least one sensor competing for access with probability 
Pcow (39). 

Orriss et al. (2002) showed that the number of sensors uniformly distributed on an infinite 
plane that hear one particular sink as the one with the strongest signal power (i.e., the number 
of sensors competing for access to such sink), is Poisson distributed with mean 


1—e r 


, 40 
Ho e) 


S 
| 
E 
un 


with us = ps Ac, being the mean number of sensors that are audible by a given sink. Such a 
result is relevant toward our goal even though it was derived on the infinite plane. In fact, 
when border effects are negligible (i.e., Ag, < 0.1A) and k is large, n can still be considered 
Poisson distributed. The only two things that change are: 


e nis upper bounded by k (i.e., the pdf is truncated) 


e the density ps is to be computed as the ratio k/ A [m ?], thus yielding ps; = ka f 


Therefore, we assume n ~ Poisson(ñ), with 


Ag. 1— g sink Ls eg |14«/À 
TE a ci olea MENS it e (41) 
l ) A Hsink I 


Finally, by taking the average in (37) explicit and neglecting border effects (see (39)), we get 


—IA,./A 1 k fe^ 
Ps = (1 =g Os ) . M X PyAc(n) ni y (42) 
n=1 ` 
where " 
zn o5—n 
M=} T (43) 


is a normalizing factor. 


7.2 Area Throughput 

The amount of samples generated by the network as response to a given query is equal to 
the number of sensors, k, that are present and active when the query is received. As a conse- 
quence, the average number of data samples-per-query generated by the network is the mean 
number of sensors, k, in the observed area. 

Now denote by G the available area throughput, that is the average number of samples gen- 
erated per unit of time, given by 


i) 
| 
A 
=> 
| 


ps: A: T [samples/sec]. (44) 
q 


From (44) we have k — GT}. 
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The average amount of samples received by the infrastructure per unit of time (area through- 
put), S, is given by: 


-Foo 
S = y S(k)- gy [samples/sec], (45) 
k=0 
where " 
S(k) = 7 Pale (46) 


gx as in (1) and Pi as in (42). 
Finally, by means of (42), (43) and (44), equation (45) may be rewritten as 


1—e-14s/A 
ga ius cR 
Tq 
AEEXaPuacQU g^ (GTy)ke OM m 
= x P (k EN 1)! 


7.3 Numerical Results 

In this section the area throughput obtained with the two modalities Beacon- and Non Beacon- 
Enabled, considering different values of D, SO, Nets, T; and different connectivity levels, is 
shown. 
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Fig. 9. S as a function of G, for the Beacon- and Non Beacon-Enabled cases, by varying SO, 
Ngrs and T}, having fixed D = 10. 


In Figure 9, S as a function of G, when varying SO, Nets and T; for D = 10, is shown. The 
input parameters that we entered give a connection probability Pcon = 0.89. It can be noted 
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that, once SO is fixed (Beacon-Enabled case), an increase of Ngrs results in an increment of 
S, since Py ac increases. Moreover, once Ngrs is fixed, there exists a value of SO maximising 
S. We can note that, a part for the case, Beacon-Enabled with GTSs allocated, an increase of 
SO results in a decrement of S. In fact, even though Pmac gets greater the query interval 
increases and the number of samples per second received by the sink decreases. On the other 
hand, when the Beacon-Enabled mode is used and GTSs are allocated, the optimum value of 
SO is 1. This is due to the fact that, having large packets, when SO — 0 too many packets are 
lost, owing to the short duration of the superframe. 

Concerning the Non Beacon-Enabled case, in both Figures it can be noted that, by decreasing 
Tj, S gets larger even though Py Ac decreases, since, once again, the MAC losses are balanced 
by larger values of fy. 
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Fig. 10. S as a function of G, in the non beacon-enabled case, for different values of D and 
Pcon, having fixed T; to the maximum delay. 


Finally, we show the effects of connectivity on the area throughput. When Pcoy is less than 
1, only a fraction of the deployed nodes has a sink in its vicinity. In particular, an average 
number, k = PcoNGT;/ I, of sensors compete for access at each sink. In Figure 10 we consider 
the non beacon-enabled case with D = 2, T; = 128 T and D = 10, T; = 136 T. When D = 10, 
T; = 136 T, for high G the area throughput tends to decay, since packet collisions dominate. 
Hence, by moving from Pcoy = 1 to Pcon = 0.89, we observe a slight improvement due to 
the fact that a smaller average number of sensors tries to connect to the same sink. Conversely, 
when D = 2, T; = 128T, S is still increasing with G, then by moving from Pcow = 1 to 
Pcon = 0.89, we just reduce the useful traffic. Furthermore, when Pcon = 0.15, the available 
area throughput is very light, so that we are working in the region where Pmac(D = 2, T; = 
128T) < PyAc(D = 10, T; = 136 T), resulting in a slightly better performance of the case with 
D = 2. Thus we conclude that the effect of lowering Pcon results in a stretch of the curves 
reported in the previous plots. 
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9. List of acronyms 

r.v. random variable 

PAN Personal Area Network 

CAP Contention Access Period 

CFP Contention Free Period 

CSMA carrier-sense multiple access 
CSMA/CA carrier-sense multiple access with collision avoidance 
GTS Guaranteed Time Slot 

ISM industrial scientific medical 

MAC medium access control 

p.d.f. probability distribution function 
PPP Poisson Point Process 

PAN personal area network 


WSN wireless sensor network 
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